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Abstract 

The pressure dependence of the Born effective charge, dielectric constant and 
zone-center LO and TO phonons have been determined for 3C-SiC by a hnear 
response method based on the hnearized augmented plane wave calculations 
within the local density approximation. The Born effective charges are found 
to increase nearly linearly with decreasing volume down to the smallest volume 
studied, V/Vq = 0.78, corresponding to a pressure of about 0.8 Mbar. This 
seems to be in contradiction with the conclusion of the turnover behavior 
recently reported by Liu and Vohra [Phys. Rev. Lett. 72, 4105 (1994)] for 6H- 
SiC. Reanalyzing their procedure to extract the pressure dependence of the 
Born effective charges, we suggest that the turnover behavior they obtained 
is due to approximations in the assumed pressure dependence of the dielectric 
constant Eoo, the use of a singular set of experimental data for the equation 
of state, and the uncertainty in measured phonon frequencies, especially at 
high pressure. 
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Typeset using REVT^ 



I. INTRODUCTION 



Recently, Liu and Vohra [jT| presented intriguing evidence regarding the pressure depen- 
dence of the Born effective charge in 6if-SiC. These authors found that the effective charge 
increased initially with increasing pressure, reaching a maximum at about 0.4 Mbar. Fur- 
ther increasing the pressure, however, resulted in a decrease in the magnitude of the effective 
charge. Because of its great potential in device applications, especially in harsh environ- 
ments, it is important to accurately characterize the properties of this unusual material. SiC 
crystallizes in hundreds of polytypes, corresponding to different stacking sequences of Si-C 
bilayers in the cubic [1 1 1] direction. Polytypes with cubic, hexagonal and rhombohedral 
symmetry are designated 3C, nH and nR respectively, where n is the number of layers in 
the repeating unit, with the 3C structure corresponding to the zincblende structure. The 
different polytypes have very similar properties, since differences in local atomic coordina- 
tion first appear in the second neighbor shell and the experimental equation of state for the 
different polytypes are very similar. Phonon dispersion along the stacking direction is 
considered universal for different prototypes and this universality has been used to map out 
the phonon frequencies in this direction from Raman scattering measurements. |^ Indeed, 
a recent self-consistent study of 3C, 2H and AH SiC by Karch et al. found only small 
differences in the calculated phonon frequencies along the stacking direction. The main 
purpose of the present work is to report a first-principles study of structural and dynamical 
properties of SiC under high pressure. In particular, the volume dependence of the Born 
effective charges and the dielectric constant were obtained. We find that the Born effective 
charge increases nearly linearly as the volume is decreased. This is in sharp contrast with 
the above-mentioned turn-over behavior reported recently for 6iJ-SiC. The likely reasons 
for this discrepancy will be discussed in light of the calculated results. 
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II. METHODOLOGY 



First principles total-energy calculations in the local-density approximation (LDA) were 
performed using the linearized augmented plane wave (LAPW) method [0] to determine 
the ground state structural properties of 3C-SiC at equilibrium and under pressure. Lat- 
tice dynamical properties, Born effective charges, and the dielectric constant were obtained 
with a linear response algorithm developed recently within the LAPW method. [§|-|T^ To 
dispense with the need to treat the chemically inert localized inner core orbitals, we employ 
a hard Kerker [|TT| type pseudopotential. The muffin-tin radii for Si and C are 1.79 and 
1.50 a.u., respectively. We used the Wigner interpolation formula for exchange-correlation 
potential. [|l^ Born effective charges and the dielectric constant were calculated using up to 



a 12 X 12 X 12 uniform Appoint mesh |T^, which corresponds to 28 special points in the ir- 
reducible Brillouin zone (BZ). At this level of convergence, the acoustic sum-rule is satisfied 
with an error of only 0.6%. In calculating the dynamical matrices, we used an 8 x 8 x 8 fc- 
point mesh. Tests with denser meshes show that phonon frequencies are converged to better 
than about 0.5%. Tests with other exchange-correlation functionals also display differences 
at this level. Phonon dispersions in the harmonic approximation were obtained in the full 
BZ as follows. First, ab initio calculations were carried out to determine the dynamical 
matrix at the irreducible phonon wavevectors corresponding to a 4x4x4 uniform mesh, 
which by symmetry gives the dynamical matrix at all mesh points. The dynamical matrix 
can then be obtained at arbitrary wavevectors by interpolation, first separating it into a 
long-range dipole-dipole term and a short-range term. 0] The former is evaluated exactly 
from the calculated Born effective charges and dielectric constant using the Ewald summa- 
tion technique. The remaining short-range part is then interpolated using real-space force 
constants, which are found through Fourier transform. 



III. RESULTS AND DISCUSSION 
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A. Structural Properties, Equation of State and Phase Transition 



LAPW total energy calculations were performed to determine the theoretical lattice 
constant Oq, bulk modulus Bq and its pressure derivative B'q for 3C-SiC by fitting to the 



Murnaghan equation of state. |15| Earlier experimental and theoretical studies suggest that 
the pressure-volume relations for different polytypes are similar, due to the rigidity of the 
nearest-neighbor coordination. The equation-of-state data for the two polytypes 3C-SiC and 
6i^f-SiC have been found to be essentially the same up to the transition pressure of 3C-SiC. 
Q Theoretical calculations |^ also yield very small difference in the bulk modulus (less than 
1%) between the 3C, 2H and AH polytypes. The calculated pressure derivative of the bulk 
modulus for these polytypes are different by only 7%. 

Results for 3C-SiC are presented in Table | and compared with other calculations and 
experiment for 3C and other polytypes. Our results agree well with other calculations given 
in the table. Slight differences may be due to the use of different forms of the exchange- 
correlation potential. There is nearly perfect agreement between the present calculations and 
those of Chang and Cohen [|16|, both using the Wigner exchange-correlation potential. Karch 
et al. used the Ceperly and Alder exchange-correlation formula [|1^ as parametrized by 



Perdew and Zunger. |T8[ Lambrecht et al. |]T9[ used an LMTO basis set and the von Barth- 
Hedin parametrization of the exchange-correlation energy. |^ As also seen in other systems 



21| , the Wigner form tends to yield slightly larger equilibrium lattice constants than other 



forms, and as expected the bulk modulus is accordingly slightly smaller. 

A number of high-pressure experiments have been carried out on various polytypes of 
SiC. While they generally show little difference between the equation of state for different 
polytypes, there are large discrepancies between the results of different experimental groups 
for the bulk modulus and its pressure derivative. Most earlier experiments report a bulk 
modulus around 2.25 Mbar. The recent work of Yoshida et al. on SiC polytypes up to 0.95 
Mbar, the highest pressurization on these materials so far, reported a relatively large bulk 
modulus Bq (2.60±0.09 Mar) and a smaller pressure derivative B'q (2.9±0.3). They ascribed 
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the discrepancy with the previous measurements to fitting to the larger pressure range in 
their experiment. Analysis of our calculated results does not support this explanation. 
Using our calculated total-energy at various volumes, we can fit to the Murnaghan equation 



of state [jT5[ for pressures up to about 0.8 Mbar. We find only small differences (both Bq and 
B'q change about 1 %) compared to fitting over a smaller pressure range. The explanation 
is also inconsistent with the experiments of Aleksandrov et al. and Goncharov et al. 
1^, which were carried out respectively to 0.425 Mbar and 0.45 Mbar but obtained smaller 



bulk modulus as compared with that of Strossner et al. ||2^, which went up to only 0.25 
Mbar. Based on the results of all existing well-converged theoretical calculations, it seems 
reasonable to suggest that the lower range of the experimental bulk modulus is more likely 
to be correct. 

We have also determined the volume dependence of the total energy for rocksalt struc- 
ture SiC. Table || compares the calculated transition parameters to other calculations and 
experiment. The agreement between the different calculations is good, predicting that the 
transformation from the zincblende phase to the rocksalt phase takes place at about 0.65 
Mbar. The volume of the zincblende phase at the phase transition is predicted to be about 
82% of its equilibrium volume at ambient pressure, and the volume reduction accompanied 
by the transition is predicted to be about 20%. The predicted transition pressure and the 
volume of the 3C-SiC phase before the transition differ considerably from the experimental 
data, but the observed volume reduction is in good agreement. This may be due to the 
fact that the experimental transition pressure was obtained for the forward transition (from 
zincblende phase to rocksalt phase), where an excess pressure beyond the equilibrium value 
appears and is included in the measurements. [H[ 
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B. Born Effective Charge, Dielectric Constant, Phonon Frequencies, and Elastic 

Constants at Equilibrium Volume 



The Born effective charge Z* and dielectric constant are calculated at small q wavevec- 
tor 0.01 [1 1 1]^, where a is the lattice constant. Table |TT1| presents our results at the ex- 
perimental equilibrium volume. The Born effective charges Z*(Si) and Z*(C) are converged 
to better than one percent when BZ integrations are performed using a 28 special /c-point 
set in the irreducible BZ, corresponding to a 12 x 12 x 12 uniform A;-point sampling. The 
acoustic sum rule requires that Z*(Si) and Z*(C) have the same magnitude and opposite 
sign. Violations arise from finite fc-point sampling. However, good results can be obtained 
with 10 special k points, especially after averaging the magnitudes of -Z'*(Si) and Z*(C). 
The averaging also reduces the effect of using a small but finite wavevector. This effect was 
checked by using a slightly larger wavevector q=0.02 [1 1 1]^, which yielded only a slightly 
larger discrepancy. Our calculations give a Born effective charge in excellent agreement with 
experiment and with the linear response calculations by Karch et al, using a plane wave 
basis set . The static dielectric constant eoo is about 7% too large, however. The tendency 
of the LDA to overestimate eoo is a well known problem that is attributed not to the density 
functional theory, but to the LDA used for the exchange-correlation potential. [pT] , p2[ 

Our calculated phonon frequencies are shown in Fig. |1|. For points along the (^] 
direction, we obtained the dynamical matrices as described in Section II, using a uniform 
4x4x4 phonon wavevector mesh. But along the 0] and ^ 0] directions, we have 
calculated the dynamical matrices at additional points, which correspond to an 8 x 8 x 8 
mesh, and obtained one-dimensional interplanar force constants to perform the interpolation. 
||33| This improves the interpolated dispersion in the acoustic region. No neutron scattering 
data for zincblende SiC is available because it is difficult to grow large enough 3C-SiC 
single crystals. The filled squares in the figure are experimental phonon frequency data 
from Raman measurements. The experimental results on the BZ boundary (X and L) are 
obtained from second-order Raman scattering . Others are first-order Raman data from 
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hexagonal and rhombohedral polytypes that have been unfolded into the larger 3C-SiC BZ. 
The excellent agreement shows that the different stacking sequences have little influence 
on the vibrational properties. As mentioned, the calculated phonon frequencies of Karch 
et al. on 2H- and 4i/-SiC directly support this conclusion. Table |V| compares our 
calculated phonon frequencies at high-symmetry /c-points with those of Karch et al. |^ and 
with experiment. The zone-center cjlo frequency is obtained from the calculated wto, Coo, 
Z* , and the volume of the primitive unit cell V , using the relation 

2 2 47re2Z*2 

-LO = -^TO + ^;^, (1) 

where /i is the reduced mass. Our calculated frequencies are systematically about 1% smaller 
than those of Ref . @] . This is due in part to the difference in the equilibrium volume in the 
two calculations. 

The elastic constants can be extracted from the calculated acoustic phonon dispersions. 
For the cubic structure, the elastic constants C44 and Cn are simply related to the velocities 
of transverse and longitudinal plane waves in [1 0] direction as follows 



..[100] = (2) 



^.[100] = ^^, (3) 

where p is the mass density of SiC. The velocities of i^t[100] and t;L[100] can be obtained 
from the slopes of the acoustic transverse and longitudinal phonon dispersions at the zone 
center in the [1 0] direction: 

^.[100] = % (4) 

where uoj is the angular frequency of the j-th branch phonon and k is the magnitude of 
the corresponding wavevector. Similarly, the elastic constant C12 can be obtained from the 
relationship |^ 

„,|nol^/^. (5) 



where vt^ is the velocity of transverse wave with atomic displacement along [1 10], which 
corresponds to the lowest acoustic branch in the S direction, i.e. [110]. The bulk modulus 
Bq is related to these elastic constants by the relation pSf 



Cii + 2Ci2 , s 

Bo = • (6) 

The directional dependence of the transverse sound velocity in 3C-SiC is related to the fact 
that the interatomic forces in 3C-SiC are highly noncentral, which implies that the elastic 
constants do not satisfy the Cauchy relation (C12 = C44). Table ^ presents the elastic 
constants and the bulk modulus as extracted using the equations above, together with other 
calculated results and experiment. Our results are in good agreement with experiment and 
slightly smaller than the calculations by Karch et al. due again in part to the use of our 
slightly larger lattice constant. The bulk modulus Bq determined from the elastic constants 
is consistent with our total energy calculations in the previous section. 

C. Pressure Dependence of Phonon Frequencies, Dielectric Constant, and Born 

Effective Charge 

The mode Griineisen parameters, defined as 

V du 

describes the volume dependence of the phonon frequencies. The calculated mode Griineisen 



parameters at the F, X and L points in the BZ are presented in Table Our results 
appear to be in good agreement with Karch et al. |^ as presented in their Fig. 9. In 
the text of that paper, they also gave the numerical value for the TA{X) mode, 0.13, 
which compares favorably with our value of 0.12. The experimental results of Olego et al. 



3^ , |36[1 and of Aleksandrov et al. ||26| are presented for comparison. In the former work, the 
Griineisen parameters were calculated from their experimental volume dependence using a 
bulk modulus of 3.219 Mbar. This value is hkely to be too large, as discussed above. The 



8 



values given in Table ^ are converted from their results, assuming the bulk modulus from 
our calculations. The results are seen to be in generally good agreement with our calculation. 
A significant discrepancy exists for the LA mode at the L-point. The experiment of Olego 
et al. gave a negative value, in contrast with our theoretical prediction and that of Karch 
et al. 10 

Fig. 1^ compares the calculated pressure dependence of the zone-center phonon frequencies 
of 3C-SiC with the experiments for both 3C-SiC and 6if-SiC. The numerical values of our 



frequencies are presented in Table [VI I|, together with the calculated Born effective charge and 
dielectric constant of 3C-SiC at different lattice constants. The second row of Table |VlT| gives 
the results at the equilibrium volume and the smallest lattice parameter corresponds to a 
pressure of about 0.8 Mbar. The calculated pressure as a function of volume was determined 
by fitting the volume dependence of the total energy to the Murnaghan equation of state 



||15|| . We have also tested fitting the Birch-Murnaghan equation of state |3^, which was 
used by Yoshida et al. [Q, and find only small differences, with Bq increasing by 1%, B'q 
increasing by less than 3%, and the pressure values changing less than 1%. As in the ambient 
pressure case mentioned above, our phonon frequencies are nearly uniformly a few percent 
smaller than those of 3C-SiC measured up to 0.25 Mbar by Olego et al. [^]. Therefore the 
experimental pressure dependence of phonon frequencies of zone-center LO and TO modes 
are fairly well reproduced in our calculations. This is consistent with the agreement of zone- 



center mode Griineisen parameters in Table VI between theory and experiment. Fig. also 



compares our calculations with Liu and Vohra's measurements for QH-'&iC With increasing 
pressure, our calculated frequencies stiffen slightly faster than their experiment, with values 
initially about 2% smaller than measurements changing to 2% larger at the end. To see 
the differences more clearly, the resulting LO-TO splitting from theory and experiments 
is presented in Fig. ^ For 3C-SiC, Olego et al.'s measured results under pressure |]3^ 



appear to increase faster than our calculations. In QH-SiC, Liu and Vohra ||T| observed a 
rapidly increasing LO-TO splitting at low pressure which saturates at the high end of the 
pressure range studied (0.9 Mbar). Our calculations compare well with the trend of the 
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pressure dependence in their measurements, but the calculated values are smaller due to the 
overestimation in the dielectric constant. The pressure dependence of this overestimation is 
unknown, of course. However, if we assume that this is not a significant factor and reduce 
the calculated £00 by 7%, the calculated LO-TO splitting increases as shown by the dotted 
lines in Fig. ^, which is in good agreement with the result of Liu and Vohra, except at the 
highest pressures. 

As mentioned, using Eq. (||), Liu and Vohra [0] deduced the pressure dependence of 
the Born effective charges from the LO-TO splitting. Their results show that the effective 
charge increases at low pressure but reaches a maximum at about 0.4 Mbar and decreases 
as the pressure is increased further. Actually, the Born effective charges in the hexagonal 
polytypes are slightly anisotropic and are different for atomic displacement in the planes 
and perpendicular to the planes. The difference is about 7% for 2H and AH polytypes. 
This anisotropy in the effective charge and other physical quantities (such as the dielectric 
constant) was neglected in the analysis of Liu and Vohra. [|1| If we consider the averaged 
effective charge, the results of Karch et al. |^ show that there is little distinction between 
the cubic and the hexagonal polytypes. In this sense, the results of our calculation for the 
cubic SiC, discussed below, are pertinent to the above-mentioned experiment. 

In order to extract the pressure dependence of the effective charge from the measured 
LO-TO splitting, two additional pieces of information are required: i) the equation of state 
to derive the volume from the measured pressure, and ii) the volume dependence of Soo- By 
contrast, our calculated Born effective charge is determined directly and is known to be very 
accurate, unlike the dielectric constant which suffers from the LDA approximation. For their 
equation of state, Liu and Vohra used the values of Yoshida et al. [|] for the bulk modulus, 
B = 2.60 Mbar, and its pressure derivative, B' = 2.9 together with the Birch-Murnaghan 
equation of state. However, this bulk modulus is 10 - 20% larger than the other measured 
values, and the pressure derivative is about 20 - 35% smaller. 

The volume dependence of Soo has apparently never been measured for SiC. Liu and 



Vohra followed Olego et al. in using the following expression, 
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^^=r (8) 
d\nV ' ^ ' 

with r = 0.6 for this range of pressures, citing earher work on both Si and C that reported 



this value. Goni et al. have reviewed the volume dependence of the refractive index of 
Ge and GaAs. Their Table II shows large differences between various measurements of the 
volume derivative of ^oo in Eq- (H)- Thus, the use of the scaling relation, Eq. (||), and the 
particular value r = 0.6 must be regarded as a weak link in the experimental analysis of the 
volume dependence of the Born effective charge. Indeed, the present calculations suggest 
that this value of the exponent scaling power is too large. In Fig. ^, we compare the scaling 
relation Eq. (||) with the calculated volume dependence of £00- The value r ~ 0.3 is found 
to give the best fit to the calculated values. As mentioned, the LDA tends to overestimate 
£005 in this case by about 7%. The volume dependence of this error is unknown, but it would 
have to be large to modify the above conclusions. 

In Fig. we compare the volume dependence of the Born effective charge Z* from 
calculations and those derived from experiment. The calculated Born effective charge for 
3C-SiC is seen to linearly increase with decreasing volume (see curve (a)), and there is 
no evidence of a relative maximum as is in the experimental results of Liu and Vohra for 
6i^f-SiC [in . The other four curves are all derived from experiment, using Liu and Vohra's 
quadratic fit to the measured LO-TO splitting for 6/f-SiC, but with different choices of 
equation of state and volume dependence of dielectric constant. In curve (6), we have used 
the same input data as Liu and Vohra [|T|: the scaling relation Eq. (§) with r = 0.6 (using 
the experimental £00 = 6.52 at ambient pressure) and the equation of state {B = 2.60 Mbar 
and B' = 2.9) from Yoshida et al. This curve exhibits the same peaked behavior as that 
shown in Fig. 3 of Liu and Vohra Curve (c), which is plotted with the same input data as 
for curve (6) except using our calculated equation of state {B = 2.20 Mbar and B' = 3.71), 
is much less peaked than curve (6), indicating that the volume dependence of Z* derived 
from experiment following Eq. (|l|) is quite sensitive to the equation of state. The other two 
curves [d) and (e) are derived using the best scaling relation Eq. (0) with r = 0.3 (again 
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using the experimental Soo = 6.52 at ambient pressure). We have used the equation of state 
from Yoshida et al. for curve (rf) and our calculated equation of state (for 3C-SiC) for curve 
(e). These two curves give larger Born effective charges than (6) and (c) since the scaling 
relation r = 0.3 yields larger dielectric constant under pressure than r = 0.6. Curve (e) is 
seen to agree quite well with the calculated Z* except for volumes smaller than about 0.85 
Vq. Examination of Liu and Vohra's measurements in Fig. which shows the scatter in the 
experimental values of the LO-TO splitting, indicates that the greatest uncertainties occur 
for the higher pressures (smaller volumes) due to the pressure induced broadness of Raman 
and ruby peaks The magnitude of this scatter ~ 4:cmr^ (or 2%) taken together with 
the questionable use of the scaling relation Eq. (|D with r = 0.6 and the large bulk modulus 
of Yoshida et al. undermines the conclusions of Liu and Vohra regarding the peaked 
behavior of Z* . Given the uncertainties in the experimental analysis, we would assert that 
our prediction of a linear increase of Z* with decreasing volume, up to pressures of 0.8 Mbar 
is not contradicted by their measurements. 



IV. CONCLUSIONS 

We have presented a detailed analysis of the volume dependence of various properties of 
3C-SiC. The calculated equilibrium lattice constant, bulk modulus and pressure derivative 
of the bulk modulus of 3C-SiC were found to agree well with experiment and other LDA 
calculations. The transition pressure from the zincblende phase to the rocksalt phase is de- 
termined to be 0.65 Mbar, in good agreement with other theoretical results but lower than 
the reported experimental value (1.00 Mbar). The volume reduction, however, is in good 
agreement. The discrepancy may be due in part to the fact that the experimental tran- 
sition pressure obtained in a diamond-anvil apparatus is for the forward transition (from 
zincblende phase to rocksalt phase), which requires an excess pressure beyond the equi- 
librium value. The equilibrium calculated Born effective charge Z* is in excellent agreement 
with experiment, but the dielectric constant is about 7% too large, as is typical in LDA 
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calculations. The phonon dispersion and the elastic constants agree with available experi- 
ment data, demonstrating that the different stacking sequences in silicon carbides have little 
inffuence on the vibrational properties. The experimental Griineisen parameters are well 
reproduced, and predictions are given for those experimentally unavailable. In 3C-SiC, we 
find that i) the dielectric constant ^oo decreases with increasing pressure, and ii) the Born 
effective charge increases monotonically with pressure, without any evidence of a relative 
maximum as reported for 6if-SiC by Liu and Vohra After analyzing the procedure used 
to extract the experimental Born effective charges Z* of 6if-SiC in Ref . , we suggest that 
the turn-over behavior of Z* reported by Ref. [jll is due to the assumptions regarding the 
volume dependence of dielectric constant ^oo, the use of a singular set of experiment data 
for Bq and Sq, and uncertainties in the measured phonon frequencies, especially at high 
pressure. Our calculations predict a linear increase of Z* with decreasing volume, up to 
pressures of 0.8 Mbar. 

Note Added. Since the submission of the manuscript, two recent publications P0|J^ 



have come to our attention. A theoretical study using molecular dynamics simulations with 
an empirical potential model EDI raised the possibility of amorphization of SiC under high 
pressure. This seems to be inconsistent with the experiment of Yoshida et ai, [^] which 
observed a direct polymorphic transition to the rocksalt structure at 1 Mbar. Thus further 



experimental and theoretical work needs to be done to see if this is a real effect. Ref. |41 



presented calculated Born effective charges of 3C-SiC under low pressure, in good agreement 
with this work 
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TABLES 

TABLE L Equilibrium state properties of SiC. The equilibrium lattice constant, bulk modulus 
and its pressure derivative are represented by oq, Bq, and B'q, respectively. 



(A) 



Bo (Mbar) 



B' 



Pressure range (Mbar) 



Theory 
Chang et al."" 
Lambrecht et al} 
Karch et al.^ 
This work 



4.361 
4.315 
4.345 
4.360 



2.12 
2.23 
2.22 
2.10 



3.7 
3.8 
3.88 
3.71 



Experiment 
3C 



3C 
3C 

6H,15R 
3C,6H 



4.360*^ 



2.24^ 

2.234^ 

2.253 

2.27± O.OS'' 
2.48± 0.09* 
2.24± 0.03^' 
2.60±0.09''' 



4.1 ± 0.10 
4.0 ± 0.3' 
4.3± 0.3^' 
2.9±0.3'' 



h 



0.425^" 
0.25* 
0.45^' 
0.95^ 



* Performed on specimens of a-SiC (hexagonal and rhombohedral phases). 
" ReL 0, * ReL 0, " ReL §, ReL H], ^ Estimated for 3C-SiC by Yean et al. | 
Schreiber et al. . 



9 R. D. Carnahan P 
^ Aleksandrov et al. |2 
* Strossner et al. p7|| . 
^ Goncharov et al. 0. 
^ Yoshida et al. M. 
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TABLE II. Transition parameters from zincblende structure to rocksalt structure for SiC. Vt 
is the volume of zincblende (3C) phase at transition, and Vq is the volume at ambient pressure. 
AV/Vq denotes the percentage of volume reduction when the transition occurs. 



Pressure (Mbar) 



Vt/Vo 



AV/Vo 



This work 
Chang et al."" 
Cheong et al} 
Christensen et al.^ 
Experiment^ 



0.65 
0.66 
0.60 
0.59 

1.00 



0.817 
0.81 
0.825 
0.84*^ 
0.757 



20.2% 
18.5% 



19%^^ 
20.3% 



Ref. fiel. 



Ref. m 



^ Ref. dg. 
Quoted from Ref. 



in which the authors derived these numbers from Fig. 12 of Ref. [29 



Ref. 
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TABLE III. Calculated Born effective charges Z* and eoo at experimental volume. 





Z*(Si) 


Z*(C) 


|Z*(Si)| + iZ*(C)| 
2 




This work 










Number of special /c-points 










2 


2.428 


-3.392 


2.910 


9.426 


10 


2.701 


-2.716 


2.708 


7.137 


28 


2.709 


-2.693 


2.701 


7.005 


Karch et al.°' 










(at LDA volume) 






2.72 


6.97 


Experiment 






2.697'' 


6.52^ 



^ Ref. [|. 
^ Ref. [g. 
^ Ref. [g2|. 
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TABLE IV. Comparison of phonon frequencies (cm ^) in 3C-SiC. 







Tlo 










Lta 


Lla 


Lto 


L-LO 


This work 


774 


945 


361 


622 


741 


807 


257 


601 


747 


817 


Karch et al.^ 


783 


956 


366 


629 


755 


829 


261 


610 


766 


838 


Experiment^ 


796 


972 


373 


640 


761 


829 


266 


610 


766 


838 



'^Reference |6[. 
^Reference [^. 
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TABLE V. Elastic constants and bulk modulus (in Mbar) of 3(7-SiC. 



C44 Cii C12 Bq 

This work 2.41 3.84 1.32 2.16 

Karch et al."- 2.53 3.90 1.34 2.19 

Lambrecht et al} 2.87 4.20 1.26 2.24 

Experiment^ 2.56 3.90 1.42 2.25 



« Ref. § 
Ref. 



The experimental data are derived by Lambrecht et al. [|T^ from the sound velocities of 
Feldman et al. ISll. 
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TABLE VI. Griineisen parameters of phonon modes in 3C-SiC. 





± ills WOi K 




LU{L ) 


i.U ( 


1 no'^ 1 1 no'' 


LO(r) 


1.02 


1.01 , 1.091 


TA ( Y\ 


n 1 9 




J \ f Y\ 


U.oZ 




iU(A j 


i.4D 


i.oU 


LO(X) 


1.16 




TA(L) 


-0.13 


-0.28° 


LA(L) 


0.90 


-0.11'' 


TO(L) 


1.31 


1.24'^ 


LO(L) 


1.15 


1.30" 



" Converted from the results of Olego et al. [^^; see text 



Aleksandrov et al. [2(: . 
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TABLE VII. Volume dependence of Born effective charge Z, dielectric constant £00 and LO-TO 



iplilliiig in 3C-S 


iC. TIk; lallice 


prU'amelcr a is in a. 11. 


and frcqnc^ncic^s 


are in cm 




a 


Z* 


^00 


wto 




<^LO — ^TO 


8.339 


2.676 


7.111 


744 


910 


166 


8.239 


2.701 


7.005 


775 


945 


171 


8.157 


2.721 


6.927 


800 


975 


175 


8.069 


2.744 


6.852 


831 


1010 


179 


8.000 


2.761 


6.800 


855 


1038 


183 


7.864 


2.793 


6.707 


905 


1094 


189 


7.710 


2.830 


6.626 


965 


1161 


196 


7.610 


2.850 


6.575 


1015 


1214 


199 
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FIGURES 

FIG. 1. Calculated phonon frequencies of 3C-SiC at the experimental lattice constant (solid 
lines). The solid square symbols are from Raman measurements (see text). 

FIG. 2. Pressure dependence of zone-center phonon frequencies of SiC. Filled circles are calcu- 
lated results for 3C-SiC. Solid lines are the measurements for 3C-SiC up to 25 GPa (0.25 Mbar) 
by Olego et al. |^^. Triangles are the experimental results of 6if-SiC reported by Liu and Vohra 
|] to 95 GPa (0.95 Mbar). 

FIG. 3. Pressure dependence of LO-TO splitting in SiC. Filled circles connected with solid 
lines are the theoretical results for 3C-SiC using the calculated dielectric constants. Filled circles 
connected with dotted lines correspond to using the adjusted calculated dielectric constants for 



3C-SiC (see text). The dashed line represents the measurement by Olego et al. [36| for 3C-SiC 
up to 25 GPa (0.25 Mbar). Triangles are the experimental results of 6-ff-SiC reported by Liu and 
Vohra Q to 95 GPa (0.95 Mbar). 

FIG. 4. Comparison of the scaling relation in Eq. (^) to the calculated values (filled circles) 
of Eoo for several values of r: (a) r = 0.3, (b) r = 0.6, and (c) r = 0.8. 

FIG. 5. Volume dependence of the Born effective charge Z* . The filled circles connected by 
solid lines in (a) are the calculated Z* for 3C-SiC. The dotted and dashed lines (6 — e) are derived 
from Liu and Vohra's LO-TO splitting of QHS\C but using different equations of state and scaling 
relations for the volume dependence of dielectric constant (see text): (6) is determined using the 
scaling relation with r = 0.6 and Yoshida et al.^s equation of state; (c) is determined using the 
scaling relation with r = 0.6 and our calculated equation of state; (d) is determined using thescaling 
relation with r = 0.3 and Yoshida et al.'s equation of state; (e) is determined using the scaling 
relation with r = 0.3 and our calculated equation of state. 
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